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Abstract 

Based on the concept of dual cones introduced by J. Opgenorth we give an 
algorithm to compute a generating system of the group of automorphisms of an 
integral lattice endowed with a hyperbolic bilinear form. 

1 Introduction 

In his famous paper Nouvelles applications des parametres continus a la theorie 
des formes quadratiques [20] . G. F. Voronoi presented an algorithm to enumerate 
(up to scaling) all perfect quadratic forms in a given dimension n. The general 
idea for that was to compute a face-to-face tesselation of a certain cone in the 
space of symmetric endomorphisms of R n based on the pyramides induced by 
the shortest vectors of a perfect quadratic form. 

Generalizing Voronoi's ideas, M. Koecher came up with the concept of self- 
dual cones or, as he called them, positivity domains [7, 8 to obtain an alternative 
to the reduction theory of quadratic forms due to H. Minkowski. 

Another slight generalization of these ideas to so called dual cones was then 
suiting for J. Opgenorth to find an algorithm to determine a generating system 
for the normalizer ^Vgl„(z)(G) of a finite unimodular group G of degree n (cf. 
P3]), which is an essential tool e.g. dealing with crystallographic space groups 
(cf. [14]). 



"The results in this paper are partially contained in the author's master's thesis |10| written 
under the supervision of Prof. Dr. Gabriele Nebe at Lehrstuhl D fur Matematik, RWTH 
Aachen University, Templergraben 64, D-52062 Aachen, Germany 

tThe author's research is supported by the DFG Graduicrtcnkolleg 1269 "Global Structures 
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The goal here is to use these methods to derive an algorithm which gives a 
generating system of the automorphism group of an integral hyperbolic lattic^]] 
(L, $) (see Notation below). 

There exists an algorithm due to E.B. Vinberg [H] to construct the max- 
imal normal subgroup of the automorphism group of a hyperbolic lattice that 
is generated by reflections. But this algorithm terminates if and only if this 
reflection group has finite index in the full automorphism group, i.e. the lattice 
L is called reflective. But reflective lattices are very rare: For example if we 
consider the lattice Z" together with the bilinear forms induced by matrices 

:= diag(-d, 1, . . . , 1) e Z" xn for d > 0, 

then it is known, that (Z n ,£f„ ) is reflective if and only if n < 19 (cf. [T9"]V 
According to the classification of reflective hyperbolic lattices of rank 3 by D. 
Allock in Q], the highest prime divisor of the discriminant of such a lattice is 
97 (cf. [U p. 24]. Thus reflective lattices can neither occur in high dimensions 
nor for high discriminants. 

To the author's knowledge so far there is no algorithm known to determine 
the automorphism group of a general hyperbolic lattice. The algorithm pre- 
sented in this paper does at least not have theoretical limitations although in 
practice it can only handle lattices of small ranks (up to 4 in general, some 
examples of rank up to 7 work as well) and moderate discriminants. 

The paper will be organized as follows: In Section we recall the ba- 
sic definitions and key results about dual cones from [T3] which give a gen- 
eral method to determine generating systems of discontinuous groups acting 
on dual cones. The application of the results in Section on hyperbolic lat- 
tices as well as a quite powerful way to shorten the calculation time is given 
in Section 0. Some examples and statistics are presented in Section 0. These 
were calculated using the computer algebra system Magma (cf. [3 ). The 
source code for the necessary MAGMA-package AutHyp . m as well as a short 
description of the included intrinsics is available via the author's homepage 
http : //www .mi . uni-koeln. de/~mmertens. 

Notation A lattice (L, <!>) always consists of two data, a free Z-module L of 
rank n with basis B — (pi,... ,b n ) and a non-degenerate, symmetric bilinear 
form $:VxV^M where V = R ®z L. The signature of the bilinear form will 
always be given as a pair (p, —q) where p denotes the number of positive and q 
the number of negative eigenvalues of the Gram matrix of $ with respect to B 
which is denoted by B $ B := ($(&,-, bj))? d=1 . By L* we denote the dual lattice 
of L and in case that L is integral, i.e. L C £,#, we write A(L) := L# /L for the 
discriminant group of L. The automorphism group of a lattice (L, <&) is defined 
as 

Aut(L) := {g E GL(V) | Lg = L and 3>(xg,yg) = $(x,y) for all x,y E L}. 

1 ln the literature often referred to as a Lorentzian lattice 
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If (L, $) is an integral lattice, we can consider Aut(L) as a subgroup of GL n (Z) 
by fixing a basis B of L: 



where A = Note that Aut(L) acts on L from the right while it acts from 

the left on the set of Gram matrices of L. 

For any ring R let R n :— R lxn the free i?-module of rank n represented as 
a row vector. By ej we denote the ith row of the n x n unit matrix /„. 

For a subset 5 of any R- module M let (S)r submodule of M generated by 
S (mostly we will omit the subscript R if there are no confusions about the base 
ring to be worried about). Similarly, if S is a subset of some group G, we denote 
by (S) the subgroup of G generated by S. 

2 Preliminaries 

In this section we briefly recall the most important definitions and results from 
the first two sections of [TB] , 

Throughout this section, let V be a real vector space of dimension n and 



be a positive definite bilinear form on V. 

Two non-empty open subsets Vf" , V^ C V are called dual cones with re- 
spect to cr, if the following properties hold: 

(DC.l) For x e Vf , y e V 2 >0 we have a(x,y) > 0. 

(DC. 2) If Vj-° denotes the topological closure of V 4 >0 in V then for any x S 
V \ Vf there is a y E V 2 -° \ {0} such that a(x, y) < 0. The same holds 
for changed roles of x and y. 

For fixed a and dual cones Vf , V 2 >0 a discrete subset D C V 2 -° \ {0} is called 
admissible if for every x in the boundary dVi° of V^ and every e > there is 
a d G D such that a(x, d) < e (cf. [T31 Definition 1.4, Lemma 1.5]). For such a 
set D and some vector x € V^ we define (cf. [T3l Definition 1.2]) 

(i) The D -minimum of x: 



Aut(L) ^ Aut z (A) :={g£ GL„(Z) | 



cr : V x V — > R 



Hd(x) :— mina(a;, d), 



(ii) The set of D-minimal vectors of a; 



Md{x) = {d G -D | a(x,d) = 



(iii) The D-Voronoi domain of a;: 
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By [T31 Lemma 1.1] this is all well-defined. Note that the set Md(x) is always 
finite, therefore T>d{x) is a pyramide in V^~° based at the origin. 

If the Z?-Voronoi domain of a vector x G V^ contains inner points, i.e. 
dim(MD(x)) = n, then we call x a D -perfect vector. The set of all D-perfect 
vectors with D- minimum 1 shall be denoted by Pd (cf. |13[ Definition 1.2]). 
As shown in P~3] Proposition 1.8], Pd is not empty because for every y G V^ 
there exists some x G Pd such that T>D{y) Q T>d(x). The proof given there is 
constructive and can also be used to compute new -D-perfect points from old 
ones. 

Let x G Pd- A vector r G Vi \ {0} with a(r,d) > for all d G M D (x) and 
the additional property that a(r, d) = for n — 1 linearly independent vectors 
d G Md{x) is called a direction of x. Every direction r of x defines a facet of 
the D-Voronoi domain T>d{x) by 

T(r) := V D {x) n {y G V 2 | <r(r, y) = 0}. 

A direction that lies in Vf^° is called blind. 

If r is a non-blind direction, then there exists a number p > such that 
2/ := x + pr G Pd- The -D-perfect vector y is called a neighbour of x in direction 
r, the vectors x and y are called contiguous. Clearly for each D-perfect point 
there are (up to scaling) only finitely many directions, hence every x G Pd has 
finitely many neighbours. 

By p~3j Theorem 1.9], the D-Voronoi domains of the D-perfect points form a 
face-to-face tesselation of the whole cone V^ , i.e. the (undirected) D-Voronoi 
graph Yd with vertex set Pd in which two vertices are adjacent if and only if the 
corresponding _D-perfect points are contiguous, is connected and locally finite. 

Now consider a group fi < GL(V) which leaves the cone Vf* invariant and 
acts properly discontinuously on V^ , i.e. for x G Vf" 

1. the orbit xfi does not have a cluster point in Vf^° and 

2. the stabilizer fl x := Stabj^ic) := {a; G fl \ xco = x} is finite. 

For oj G ri the adjoint element of lu is the unique ui ad G GL(V) fulfilling 
o~(xuj,y) = a(x,yu> ad ) for all x,y G V. The adjoint group of fi is defined by 

n ad ._ ^ad I u g CI}. 

By [131 Lemma 2.1], the group ft acts on the Z?-Voronoi graph F^ in case 
that D is admissible. 

There is a whole theory built around groups acting on graphs, often called 
Bass-Serre theory. We refer the reader to the monographs [S] by W. Dicks and 
[17] by J. -P. Serre for the original statement of the Bass-Serre theorem, which 
implies the following theorem if translated into this context. 

Theorem 2.1. (XE3 Theorem 2.2]) 

Let f2 < GL(V) act properly discontinuously on V^ and let D C \ {0} 
discrete, admissible and invariant under the action of Sl ad . Furthermore, let X 
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denote a transversal of Po/fi such that the subgraph Trj(X) ofTrj generated by 
X is connected and put T as a maximal spanning tree ofTrj(X). The vertices 
ofTjj which are adjacent to T are collected in the set S(T). Finally, choose for 
y G <5(T) one ui y G Q such that yujy 1 G X . Then it holds that 

fl = (uj y ,Stab n {x)\y G S(T), x G X). 

In particular, is finitely generated, if the residue class graph Td/Q is finite, 
i.e. there are only finitely many D-perfect points up to action of Q. 

In case that the residue class graph is finite, the above theorem immediately 
gives a template algorithm to compute a generating system for Q: 

Algorithm 2.2. 

Input: Vi°, V^ dual cones w.r.t. a, D C \ {0} discrete, admissible and 
Q-invariant. 

Output: a transversal for Pd/Q, generators for D,. 

1. Find a D-perfect point x\ and initalize L\ = {xi}, Li := 0, S := 0. 

2. If L\ = 0, then return L2, S; else choose x G L±. 

3. Compute a generating system S x for Q x and put S :— S U S x . 

4- Compute the set R of neighbours of x and a transversal R' of R/VL X . 

5. For y G R' decide whether there is a z G £1 U L2 and an uj G f2 with 
yuj = z. 

If not, put L\ := L\ U {y}. 

If Z G h\ then put S := S U {uj}. 

6. Put L 2 := L 2 U {x}, Li := L x \ {x}. Go to step 2. 

One of the crucial problems using the above method is to prove the finiteness 
of the residue class graph Tp/Sl. For Opgenorth's normalizer algorithm this 
was proven in [5] and |12j . In general there is the following theorem which is a 
generalization of [5J Satz 6] . 

Theorem 2.3. (Koecher, 1960) 

Let L a full lattice in V and D C PI L \ {0} admissible. Assume that 
Lfl ad — L and DVl ad — D. Then the following are equivalent: 

(1) There are only finitely many SI- equivalence classes of D-perfect points with 
D -minimum 1. 

(2) There exists a set F C V^ fulfilling that for every y G V^ there is an 
ui G ft such that xui G F, such that there is a yo G V^ with d — yo G 

for all deFflL. 
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Koecher only proves the direction (1) => (2) of the equivalence for self-dual 
cones, but his proof works with some straight forward modifications for dual 
cones as well. The other direction follows immediately from an assertion (cf. [HI 
p. 400]) of which Koecher's proof c an a lso be easily adapted. 

The condition (2) in Theorem |2.3| is fulfilled for example if there is a finite 
set M C D such that the pyramide P(M) = {J2deM a dd\ad > 0} contains a 
fundamental domain for the action of fi od (cf. [131 Proposition 2.5]). 



3 Hyperbolic Lattices 

In this sec tion we are going to give detailed and explicit methods how to apply 



Algorithm |2.2| to compute generators for the automorphism group of a hyper- 
bolic lattice. For this purpose, we first have to find a pair of dual cones with 
respect to a bilinear form. 

Throughout this section, we consider integral lattices (L, $) where $ is a 
bilinear form of signature (n — 1, — 1). We also assume that there is a fixed basis 
B of L and we set A := s$ B G Z" x ™. The set of real (integral) hyperbolic 
matrices is denoted by Rj£j* (Z£*™). 

With respect to this basis, we identify V := M ®z L = R™. Therefore we will 
not always distinguish strictly between the lattice [L, $) and its Gram matrix 
A. 



3.1 Hyperbolic Lattices and Dual Cones 
Lemma 3.1. 

Consider the standard hyperbolic matrix 

H:=diag(-l,l,..A)eZ^. 

Then the set 

C := {(xi, . . ,,x n ) e K™ | xHx tr < and x x > 0} (1) 

is a selfdual cone with respect to the standard scalar product. 

Proof. Obviously, C is an open, nonempty set. 

To prove the property (DC.l), let x = (xi, . . . , x n ), y = (yi,...,y n ) £ C. 
Without loss of generality, we may assume that x\ — yi, since we can rescale y 
by a positive real number. We calculate 

n 

xy tr = x\yi + 2J %iVi 

i=2 

1 ™ 1 



G 



1 / n n n \ 

\ (4 -E^) + (vi - E »?) + E(«* - y*) 2 

\ i=2 i=2 i=2 / 



(- yj ffy tr ) + E(^ - 



> 



v >s y 

hence (DC.l) is shown. 

Regarding (DC. 2) let x = (xi, . . . , x n ) g C. If xi < 0, the j/ = (1,0, . . . ,0) € 
C fulfills xy = X\ < 0. Thus suppose iei > 0. By definition of C we know that 
— x\ + s > for s :— J27=2 x h nence especially s > 0. As above we may therefore 
assume s = 1. This yields < x\ < 1 . Now define y — (1, — x 2 , . . . , — a;„). It 
holds yHy tr — and y\ = 1 > 0, thus y £ C and xj/* 1, = xi — 1 < as desired. □ 

We remark here that for dual cones Vf* ,"^ with respect to a bilinear form 
g the sets V^°g and V^ 1 , (g ad )~ 1 are dual cones with respect to a for every 
seGL(V). 

Now Sylvester's Law of Inertia states that for any A £ ^-hyp there is a 
T £ GL„(R) such that TAT fr = H. Thus we obtain that for any such T 

Vf := CT- 1 = {x£ R n | xAx tr < and xy[ r > 0} (2) 
V 2 >0 := CT tr = {x £ W 1 | xA-V r < and xyf > 0}, (3) 

where y\ = eiT -1 and 2/2 = e\T tr , are dual cones with respect to the standard 
scalar product. 

Remark 3.2. 

Note that the dual cones Vi° and V ^ i n ([2]) and ([3]) depend on the choice of T, 
but the second output of Alaorithm \2. 6 A being a generating system for the group 
SI does not depend on the choice of the dual cones, as long as Q acts on them. 

In fact il = Autz(A) does not act on Vf* because one has —I n £ Autz(A) 
for all symmetric matrices A but certainly V("° ^ Vf* (— J n ). But Autz(^4) 
acts on Vi° U (— Vf* ) and thus we can restrict ourselves to the stabilizer of 
Vi° in Autz(.A) which from now on will be denoted by f2. We shall call fi 
the reduced automorphism group of A or L respectively. Obviously we have 
Aut z (A) = fi x (-/„). 

Lemma 3.3. 

The group fl acts properly dis continuously on V^ . 

Proof. Suppose that there is a cluster point a;* of the orbit xSl of some x £ 
V^ . Then there is a sequence (uik)k£N in fi such that xuou tends to x* as k 
tends to 00. Obviously we have xAx tr — (xLUk)A(xu>k) tr for all k, hence also 
x*Ax* r — xAx tr . Since the real automorphism group modulo — J n denoted 
acts transitively on the set {x £ Vf* | xAx tr = c} for any given number 
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c > 0, there must be an u>* e Or such that xlo* = x*. By choosing the right 
representatives modulo Stabn H (x*) we may assume that uik — > lu*, but because 
ft is discrete this means that the sequence (u>k)k£N and hence the sequence 
(a:wfe)fcGN becomes constant at some point. 

As will be shown in Lemma 1 3 . 1 2k the stabilizer of any x S V^ is isomorphic 
to a subgroup of the automorphism group of a positive definite lattice, thus it 
must be finite. 

Therefore the action of O on Vf* is properly discontinuous. □ 



3.2 Admissibility 

For the rest of this paper let 

D:=rnV 2 -°\{0}. (4) 

This set is obviously discrete and ri ad -invariant. In this subsection we want to 
prove that D is admissible, which is less obvious. 

Remark 3.4. 

For every e > and x G K n the set 

P e (x) := {y € R n | < xy tr < e} 
contains a point ^ with integral coordinates. 
Proof. For given e > and x S K™ consider 

X := {y e E" | -e < x?; tr < e}. 

This set is clearly centrally symmetric and convex and has infinite volume. 
Thus by Minkowski's Convex Body Theorem ([11] Satz 4.4]) we know that 
inZ"/ {0}. For £(E XnZ n \ {0} we either have x£ tr > and thus I G P £ (x) 
or x£ tr < and thus -I E P e (x). □ 

In addition we will need a famous result by Lovasz. To understand it we need 
the following definition (cf. [2]). 

Definition 3.5. 

(i) Let S C M" convex. We call S a maximal lattice-free convex set, for short 
MLFC-set, if it holds that 

(a) the relative interior of S with respect to (S)w does not contain points 
with integral coordinates, 

(b) S is maximal (with respect to inclusion) among all convex subsets 
T C R™ fulfilling property (a). 

(ii) Let IA an affine subspace ofW 1 . If the integral points inlA generate U as 
an affine space, i.e. 

(u nz n ) o/ f =u, 

then we call U a rational subspace ofW 1 . If not, IA is called irrational. 
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Theorem 3.6. (Lovdsz, 1989) 

Let S C W . It holds that S is an MCLF-set if and only if one of the following 
conditions hold: 

1. S is a convex polyhedron of the form S = P + C, where P is a polytope 
and C is a rational vector space, whereby it holds that dim(P) + dim£ = n 
and S° [~l Z" = and every surface of S contains an integral point in its 
relative interior. 

2. S is an irrational affine hyperplane of M n of dimension n — 1. 
A proof of this can be found in [2] . 

Remark 3.7. 

Let S = P + £ an MCLF-set as in Theorem \3.b\ U\ Then P must be a bounded 
polytope since the surfaces of S contain an integral point in their relative inte- 
rior. By adding suiting integral points if necessary one gets a basis of a full-rank 
sublattice of Z™ which has a fundamental polytope F. The projection of F onto 
(P) must be contained properly in P if P is unbounded which yields a contra- 
diction to S being lattice-free. 

We can now prove the following: 

Proposition 3.8. 

The set D from equation Q is admissible. 

Proof. Let x £ SVf . By [H Lemma 1.1] there is a y e V 2 -° \ {0} such that 
xy tr = 0. Now let £ > and define 

M £ := {v e V 2 -° \ {0} | xv tr < e}. 

We have to show that M e n D ^ 0. 

Assume on the contrary that this were false. Then, since M e is convex, there 
must be an MCLF-set S C 1™ with M e C S . O bviously, S cannot be an affine 
hyperplane, hence S = P + C as in Theorem E3HJ Since P is bounded, it can't 
hold that y € P, thus y = p + 1 for some p £ P, t G £ \ {0}. For reasons of 
dimension, p and £ are in fact unique. Now let A > 0. Thus Xy = Xp + XI, but 
since P is bounded and Ay £ M £ C S = P + £, this cannot hold for big values 
of A ifp^O. Thus y£C. 

Now consider for < e' < e the set 

T e ,(x) := {v£R n | xv tr =e'}. 
Then for each such e' we can find some v £ > £ Pi T e i{x), hence we can write 

T e i(x) = V e < + T Q (x) 

and thus 

Pe= |J (V e >+T (x)). 
0<e'<e 
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For fixed < e' < e let v e i = p' + £' where peP and £ G C. Now for every 
y G Tq(x), there is some Ao > such that 

iV + Ay G V 2 -° for all A > A . 
We can decompose v s i + Xoy = p + £ with p G P ', £ G £, thus we have 

Aoi/=(p-jO + (*-0- 
For A > Ao this implies that 

v e , + Xy = p' + A(p -p>) +£' + _ £') 
Ao Ao 
, ' » ' 

which can only hold if p — p' because P is bounded. But this means that y G C 
and hence it follows v e < + T (x) C P + £ for all s', and therefore P e (a;) C 5 



which is a contradiction to Remark 3.4. □ 



3.3 D-Minimal Vectors and Automorphisms 

Now that we have found the dual cones and an admissible set D, we give explicit 
methods to calculate the D-minimal vectors of a point x G V^ as well as its 
stabilizer and connecting elements. 

As stated before, we consider A G Z^ p n a given hyperbolic matrix with a 
fixed pair of dual cones V^ , V^ and cone test vectors 2/1,2/2 as in equations 
^ and ^ and D as in equation ^ . 

Remark 3.9. 

(i) For any x G V>° , it holds that —xA G V^°. 

(ii) Let x G Vf n Z". Then we have p D (x) < -xAx tr =: N(x). 

Proof, (i) is clear for A = H with H = diag(—l, 1, . . . , 1) and in general follows 
from Sylvester's Law of Inertia. 

(ii) follows from (i) because — xA G D and thus hd{x) < x(—xA) tr — 
-xAx tr . □ 

Lemma 3.10. 

Let x G Vf* (~l Z n and A > and define the affine hyperplane 

H\[x) := {-XxA + y e R" | xy tr = 0}. 
If A is minimal such that the finite set 

M\ = {de H x (x) n IT I (d + XxA)A- 1 (d + XxA) tr < X 2 xAx tr and y 2 d tr > 0} 
is nonempty, it holds that 

M\ = M D (x). 
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Proof. First of all we note that there actually is a minimal A such that M\ ^ 0. 
For that let d = — XxA + y G H\{x) fl Z n . Since by assumption x G Z™ and 
A G Z£*« it holds that 

a;d tr = ~XxAx^ + xy tr = XN(x) G Z, 

thus it follows that for M\ to be nonempty it is necessary that A G N \ x ^ - 
We calculate 

dA~ l d tr = ~X 2 N(x) + yA~ x y tr , 

thus minimizing A means minimizing xd tr . 

For d to belong to it is further necessary that 

(d + XxA)A- 1 (d + XxA) tr < X 2 N(x), 

and thus we get that M\ C D, and hence M\ — Md(x) for the minimal A with 
M\ ^ 0. □ 

Remark 3.11. 

There are algorithms to calculate short vectors in positive definite lattices, see 
for instance '16, pp. 188-190]. We can use them to calculate the set M\ for 
^ = jvTxJ' n\x) ' • • • un ^ -^A * s ""i empty in the following way: 

The matrix A -1 induces a positive definite scalar product on the subspace 
Hq(x) ofW 1 . We can thus calculate all y G Hq(x) such that d = — XxA+y G Z n , 
which especially means that y is contained in the lattice -jy^jZ" DHq(x) endowed 
with the bilinear form induced by A^ 1 , and 

yA- l y tr < X 2 N(x). 

Since —xA G Hi(x) is an integral point, we have to repeat this at most N(x) 
times. 



Lemma 3.12. 

(i) Since A is an integral matrix, all D -perfect vectors have integer entries 
(up to scaling). Thus we can restrict ourselves on integral vectors x. 

(ii) Let x an integral D -perfect point and g G S x := Stabn(ic). Then g induces 
an automorphism of the positive definite lattice L(x) = Hq(xA) D I/ 1 with 
scalar product induced by A. 

Proof, (i) is clear. 

Regarding (ii), we have that x G V^ and thus xAx tr < 0. In addition it is 
perpendicular with respect to the bilinear form induced by A to L(x), hence this 
bilinear form must be positive definite on the orthogonal complement Ho(xA) 
of x. That g induces an automorphism of L(x) is easily calculated. □ 
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Remark 3.13. 
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(i) There are quite powerful methods to calculate the automorphism group of 
a positive definite lattice, see for instance 115}/ . Thus we can calculate the 
automorphisms of L(x) and continue them such that they act trivially on 
x. This can be done by adding x to a lattice basis of L{x) to obtain a basis 
of M™ . With respect to this basis, the automorphisms have the form 



g e Aut(£(x)). 



Changing the basis to B chosen in the very beginning yields a subgroup 
G x of rational automorphisms of A. The elements in G x that fix Z" as a 
point set are exactly the elements in S x . 

(ii) A slight modification of this also gives a method to determine the con- 
necting elements in Theorem l-2.il . For D-perfect vectors 6 V^ with 
N(x) — N(y) we calculate integral isometries of the lattices L(x) — ¥ L(y). 
If there are no isometries, then x and y are inequivalent under £1, if there 
are then continuing them onto all of R" such that x <— > y yields again 
some rational automorphisms of A, of which it is to determine, whether it 
is integral or not. If there is an integ ral one, this yields such a connecting 



element denoted cu y in Theorem \2.i . If not, then again, x and y are not 



equivalent under fl. 

Note that this last case does in fact occur. 



3.4 Finiteness of the Residue Class Graph 

Up until now, we have given explicit methods to solve all the computati onal 
tasks in the algorithm in section |2j But it is not clear at all, that Algorithm 1 2. 2 
terminates for any integral hyperbolic lattice, i.e. for the reduced automorphism 
group of any such lattice there is only a finite number of inequivalent perfect 
points. In this subsection we will prove this using similar methods as in |18j . 
where C.L. Siegel has proven that the automorphism group of any definite or 
indefinite lattice is finitely generated. 

We recall some of Siegel's results. For that let S G Z" x ™ be an integral 
symmetric matrix of signature (to, — (n — to)) and f2 its reduced automorphism 
group. Consider the equation 



HS^H = S 



(5) 



for positive definite matrices H e IR"g" 
rational manifold f) of dimension m(n — 



H = 2Z - S, Z 



T -l 




y T -l 





, The solutions of equation ([5]) form a 
to) which can be parametrized by 



T = (I m | Y tr )S~ 



Y 



(G) 
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where where Y £ BS n m ) xm is chosen such that T is positive definite (see [TBI 
Equation (35)]). Note that fi acts from the left on Sj by 

(g,H)^gHg tr 

and from the right on the corresponding set of Y by 

{Y,g)^g tr (^h) 

where we have to renormalize the last term (by multiplication with a matrix 
Q G GL m (R)) such that the upper part becomes again I m . Thus we get a 
rational injection 

0:i3^M ( "- m)x "\ H^Y (7) 

which is fi-equivariant. 

Using reduction theory of quadratic forms due to H. Minkowski and results 
by C. Hermite, Siegel constructs a fundamental domain F of O in Sj. 

Theorem 3.14. (Siegel, 1940) 

Let S € Z nx " be an integral symmetric matrix of signature (m,—(n — trij). 
There is a finite number of hyperplanes in the subspace of K™y X m spanned by Sj 
such that ¥ is an intersection of the corresponding halfspaces and the manifold 
Sj of all matrices H fulfilling equation ([5])- In addition, ¥ has only finitely 
many neighbours, that means there are only finitely many elements ui £ f2 with 
¥H¥uj ^ 0. 

Proof. [IS, Satz 10] □ 
Equipped with these results we can now prove the main result of this subsection. 



Theorem 3.15. Let A 6 Z^ x p ra be an integral hyperbolic matrix, Q its reduced 
automorphism group and V>° and D as in equations (El) -(Til). Then the residue 
class graph Fo/fi is finite. In particular, Alaorithm \2.3i terminates. 

Proof. Obviously, we can replace A by —A without changing tt. Especially 
it holds that x(— A)~ 1 x tr > for every x £ V^ . The signature of —A is 
(1, 1 — n), thus the manifold Sj as described in equation © is homeomorphic via 
4> in equation ([7]) to the manifold V^ >0 /K>o as we can rescale each cone vector 
such that its first coordinate is 1. Since (f> is f2-equivariant, the fundamental 
domain F is m apped to a fundamental domain F C of the action of Q. 



Theorem 13.141 tells us that F and therefore F are bounded by finitely many 



hyperplanes in R n . But then Theorem 12.31 implies that there are only finitely 



many D-perfect points modulo Q, hence the assertion follows. □ 



Remark 3.16. Since Theorem \3.lA is not restricted to a certain signature of 



the matrix S , it is very likely that the method described here will yield a more 
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general procedure to determine the automorphism group of any indefinite lattice. 
In that case one would have to define the dual cones as subsets of R" xm in 
an analogous way, replacing the symbols < and > by "negative" and "positive 
definite" respectively. 

3.5 The Watson-Process 

In this subsection we are going to recall the Watson process, that enables us to 
reduce the discriminant of a quadratic form without increasing its class number. 
See e.g. [3T] and [5] for details. For that purpose, consider for a moment (L, $) 
to be an integral lattice with a symmetric, non-degenerate bilinear form <fr of 
arbitrary signature. 

Definition 3.17. 

Let p a prime number such that A(L) contains an element of order p 2 . The 
integral lattice 

FilUL) := ( P L*nL,-<f>) 
P 

is called a p-filling of L. 

The p-filling has the following effect on the genus symbol of L in the sense of 
[U Chapter 15]: For the p-filling to be defined, the highest level of a p-adic 
Jordan constituent of L must be at least 2. The p-filling reduces this level by 2. 
Repeating this until it is no longer possible yields the so called Watson lattice of 
L, denoted by Watson(L). The mapping L Watson(L) is called the Watson 
process. 

Remark 3.18. 

For an integral lattice L, the discriminant group A(Watson(L)) has squarefree 
exponent by construction and, which is even more important for computational 
issues, the discriminant |A(L)| is decreased by the Watson process by a factor 
p 2 for every p-filling during the Watson process. 

Remark 3.19. 

Since p-filling (and therefore the Watson process) cannot be expressed as an 
action o/GL„(Z), the automorphism groups of L and Fill p (L) will not be iso- 
morphic. But up to isomorphism one always has 

Aut(L) < Aut(Fill p (L)) and [Aut (Fill p (L)) : Aut(L)] < oo. 

Hence one can determine Aut(L) again via an orbit- stabilizer computation in 
Aut(Fill p (L)) or Aut(Watson(L)) respectively. 

Although we cannot prove it, so far all examples suggest that first calculating 
Autz(Watson(A) j^] via the algorithm and then finding Autz(^4) as a co-finite 
subgroup in it via an orbit-stabilizer calculation makes the computation very 
much faster, provided that the Watson process has any effect. 

2 We use the same notation for the Gram matrices of the lattices involved 
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Hcuristically, this could be explained by two observations: Since the auto- 
morphism group of Watson(L) properly contains the one of L up to isomorphism, 
thus the corresponding equivalence relation on the set of D-perfect points is in 
a way coarser for Watson(L) than for L. Thus one can expect that modulo 
Aut(Watson(L)) there will be fewer classes of -D-perfect points than modulo 
Aut(£). 

Furthermore, the computation time to determine Z?-minimal vectors of x S 
Vf n Z" depends among others on the discriminant of the lattice ( N ^ Z" n 

Hq{x),A~ 1 ) and the number N(x). Both numbers tend to be smaller for 
Watson(L) than for L. Since in most cases one has to calculate Z?-minimal 
vectors for many different vectors x, an acceleration of the computation of D- 
minimal vectors is quite valuable. 



The strength of this slight modification is illustrated in Example 14.2 



4 Examples 

All computations were executed on a Quad-Core AMD Opteron(tm) Processor 
8356 running at 1150 GHz. The used Magma version is V2.19-1. 

Throughout this section we consider the lattice (Z n ,$) where the Gram 
matrix of $ with respect to the standard basis is given by A. 

Example 4.1. 

Our first example we discuss the steps of the algorithm in detail. We take here 




A=\-3 14 8 1, A^ 1 = — I -25 12 -11 



-90 -25 10 
-25 12 -11 
10 -11 23 

For this matrix the algorithm finds 9 inequivalent D-perfect points. Rescaled 
such that their entries become integral, they are represented by 




XI = {1 


0) 


Xi 


= (2 1 


-1) 


X3 


= (2 1 


3) 


x 4 = (9 


-2) 


x 5 


= (5 3 


-3) 


x e 


= (12 5 


-7) 


x 7 = (3 


2 -1) 


Xg 


= (14 9 


-2) 


Xg 


= (21 8 


-12 



where these points each have got (in the listed order) 8, 4, 6, 8, 4, 3, 4, 3 and 6 
neighbours. The stabilizers of xi, x§, Xq, X7, Xs are trivial and all the other 
stabilizers are cyclic of order 2 given by 

/100\ /950 
Stab(zi) =( 6 -1 1) Stab(x 3 ) =( -16 -9 

\2 -1/ \-12 -6 -1, 

/125 -28 \ /1385 528 -792 \ 

Stab(ar 4 ) =( 774 -1 -172 ) Stabfag) =( 1974 751 -1128 ) 
\558 -125/ \3738 1424 -2137/ 
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In addition, we find 16 connecting elements 



c i,i 



C4,4 



C4,8 = 



'5,5 



C6,6 = 



C7,7 



cg,9 = 



561 



-5 


10 


-15 


32 


8 


-15 


-28 





84 



3920 -195 -588 
,2440 -122 365/ 

1005 5 -232N 
-1784 -9 412 
,-1048 -6 243 

51 35 -30\ 
-40 -27 24 
40 28 -23 y 

'139 63 -84' 
120 53 -72 
,320 144 -193, 

45 33 -22^ 
-56 -41 28 
8 6-3 

9 5 ON 
70 39 
,30 16 1, 

'1575 603 -902 
2506 961 -1436 
,4422 1694 -2533; 



125 -28N 
n.j = i -24 1 4 
-308 69 

'145 9 -44 
| 966 59 -292 
v 642 40 -195; 

'289 180 -180^ 
. -, = | -96 -59 60 

368 230 -229, 

109 57 -68N 
. -.,8 = | -186 -97 116 
-82 -42 51 

591 432 -224\ 
r 7.-, = | 880 643 -332 
,1620 1184 -613; 

75 51 -16> 
r„.r, | 26 17 -4 
,142 96 -29, 

231 96 -136> 
r,, 7 = | -374 -155 220 
-118 -48 69 

'1385 528 -792' 
oui I 1974 751 -1128 
,3738 1424 -2137; 



TTie notation shall be understood that via the connecting element Ci_j the D- 
perfect points X{ and Xj are adjacent in the residue class graph, in other words 
Xi and XjC^j are contiguous. The precise contiguity relations between the points 
are represented in the residue class graph Yjj/Q,. The straight lines represent 
direct contiguity, the curved ones contiguity by the connecting elements which 
are labeled on the edges. 



Calculating this information takes about 4.4 seconds. 

In the first example it does not matter whether one applies the Watson process 
since it does not have any effect on the matrix and runs in virtually no time. 
The following example shall illustrate the situation when the Watson process 
does change things. 

Example 4.2. 

In this example we want to demonstrate how effective the usage of the Watson 
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C6.6 I • 




c 9,9 



X4 J C^^,C 4 4 



Figure 1: Residue class graph Td/Q 



process can be. For that purpose, we look at the matrix 

1 




A' 



450 



41 -37 1 
-37 -16 43 
1 43 11, 



We have del A = -2 2 • 3 2 • 5 2 and the genus symbol (in the notation of J^j, 
neglecting the additional parameters of the 2-adic genus symbol) of A is given 
by 

(l-2 2 )(l- 1 -9)(l- 2 -25- 1 ), 
thus the Watson process does have an effect and yields 

Watson(A) = 




with det Watson(^) = -2 2 . 

The direct calculation of the automorphism group of A takes about 71 seconds 
while using the Watson process gives the result in next to no time. For A, the 
algorithm finds 16 inequivalent D-perfect points, for Watson(A) there is only 
one. This and the reduction of the determinant by a factor 15 2 explains this 
Terence. 



17 



4.1 Statistics 

Since we can not estimate the running time for our algorithm, we would like 
to give some statistics. For this purpose, we randomly chose 1000 integral 
hyperbolic 3x3 matrices with entries between —25 and 25 and determinant 
> —200 and measured the time (in seconds) needed to compute every single 
automorphism group, once applying the Watson process (row Watson I), once 
not (row no Watson I) , where for both tests we used the same list of matrices to 
allow direct comparison. To emphasize the strong effect of the Watson process 
we did the same thing with 1000 integral hyperbolic matrices, for which the 
Watson process really changes the matrix (rows no Watson II and Watson II). 
We also compared two runs where we required the tested matrices to be isotropic 
and anisotropic respectively (whereas we applied the Watson process in both 
cases). Table [l] gives the statistics. The value in brackets behind the average 
values is the standard deviation of the data. 





Time (sec.) 


min./max. Time (sec.) 


# Points 


min./max. # Points 


no Watson I 


95.45 (212.57) 


0.01 / 1263.99 


9.7 (9.0) 


1 / 44 


Watson I 


68.86 (195.11) 


0.01 / 1240.99 


7.6 (8.8) 


1 / 44 


no Watson II 


70.08 (178.92) 


0.02 / 1623.69 


7.4 (6.7) 


1 / 34 


Watson II 


2.96 (26.64) 


0.01 / 481.62 


2.7 (2.3) 


1 / 15 


Isotropic 


98.18 (223.73) 


0.01 / 1285.13 


8.9 (10.3) 


1 / 44 


Anisotropic 


2.43 (3.67) 


0.01 / 37.36 


5.1 (3.8) 


1 / 18 



Table 1: Computational statistics 



Apparently, automorphism groups of anisotropic lattices are - in a way - far 
easier to compute than the ones of isotropic lattices. This phenomenon occured 
while testing the implementation of our algorithm, but we cannot really explain 
it. However, this also shows the significance of the usage of the Watson process, 
because, heuristically, the Watson process affects about 50% of all isotropic 
forms of rank 3, while it only affects about 30% of all anisotropic forms and the 
computation for anisotropic forms is much faster anyway. 
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